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1. Introduction 

Recently there has been a considerable resurgence of interest in charmonium physics, moti- 
vated by the observation of a number of puzzling resonances at B factories, and by the advent of 
CLEO-c as well as by the BES upgrade. The QCD spectrum is potentially very rich, admitting 
states with high angular momenta, glueballs, four-quark states, and hybrid mesons with manifest 
gluonic degrees of freedom, and therefore numerous ideas have been advanced ^ to interpret these 
newly observed states. To resolve these ideas, a comprehensive, first-principles calculation of the 
charmonium spectrum using lattice QCD is therefore crucial. Moreover, a detailed lattice study of 
hybrid mesons in the charmonium region will help us to develop techniques that can be used in the 
light-quark sector, whose exploration is the goal of the GlueX experiment as part of the 12 GeV 
upgrade at Jefferson Laboratory. 

To date, the lattice investigations of the spectrum of exotics and states with higher angular 
momentum are quite limited. In the light-quark sector, though several calculations have been per- 
formed to study the exotic 1^+ a firm conclusion from lattice QCD about the existence of this 
state at low pion mass region is yet to be drawn. In the heavy quark sector only a handful of cal- 
culations have been performed to study hybrids and higher spin states [^, The main difficulty 
encountered when studying these states is the poor overlap of these states with the chosen inter- 
polating operators. Since it is not possible to construct an exotic-meson interpolating field from 
local quark bilinears alone, one needs to include link variables directly in the interpolating field, 
introducing additional statistical noise. Similarly, states with higher angular momentum require 
operators constructed from displaced fields, likewise introducing statistical noise. To efficiently 
tackle such states requires the use of appropriate smearing of the quark and gluon fields, and a 
large statistical ensemble. 

Besides the spectrum, it is only recently that there has been a lattice study of transition form 
factors between charmonium states [|5|]; in two exploratory studies it has been demonstrated that 
radiative decays as well as two-photon-decay form factors ^ can be studied on the lattice. We 
intend to extend those studies from form factors of conventional, low-lying charmonium states, 
to hybrids and other excited states. A subsequent study in the light quark sector will provide 
much needed information on the photocouplings to the light hybrid mesons relevant for the GlueX 
experiment. 

In these proceedings we report the preliminary results of a comprehensive study of exotics 
and higher spin mesons in the charmonium region, in the quenched approximation to QCD. We 
employ as a basis of operators those introduced in Ref. [||], together with additional operators 
missing in their construction. In forthcoming publications, we will report in more detail the results 
on the masses of ground and excited states, as well as on transition form factors between various 
states, including hybrids. These results in the charm sector will facilitate our similar study in the 
light-quark sector, using dynamical fermions. 

2. Numerical details 

2.1 Anisotropic clover action 

The relatively large value of the charm-quark mass, ~ 1.1 — 1.5 GeV, poses challenges to 
lattice QCD at the currently employed inverse lattice spacings ofa-^ ~ 2-3 GeV, since nica is then 
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large introducing substantial discretization uncertainties. Whilst a non-relativistic action might be 
appropriate for the b quark, in the charmonium sector it is more accurate to use a fully relativistic 
formulation. We employ an anisotropic lattice, with finer lattice spacing in the temporal than spatial 
directions, using a Wilson gauge action and clover fermion action. The use of an anisotropic lattice 
introduces additional parameters beyond those of the isotropic formulation, which must be tuned 
to yield the required anisotropy. We will now outline this tuning. 

In the quenched approximation to QCD, the anisotropy in the gauge sector can be tuned before, 
and independently of, the anisotropy in the fermion sector. We follow the approach of Ref. |^], 
tuning the action to yield a renormalized anisotropy = a/ao = 3, where a and ao are the spatial 
and temporal lattice spacings respectively. 

We write the fermion action in the form ^ yQY^ where 



1 V- / 1 

Q = mo + VoVo7o - -roOoAo + ^ i^W tYk - -^at^k 

^ (OoY,<^okFok + (oY,<^kiFki 

k k<l 



(2.1) 
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Detail about this action and its implementation is given in Ref. [^. Classical values of action 
parameters are given by ^ 

(2.2) 

(2.3) 
(2.4) 
(2.5) 

Following Ref. [^], classical values of the clover coefficients, Eqs. 2.4 and 2.5, are improved by 
tadpole factors, Us and Ut, for spatial and temporal links, respectively, to yield 

Cs — > Cs/u] , (2.6) 
Ct — > Ct/{utu]) . {1.1) 

The charm mass is determined by tuning the bare quark mass mo non-perturbatively such that 
the spin average of the lowest S-wave mesons (IS) coincides with its experimental value, i.e. : 
(3my/vj< + m,^J/4 = 3.067 GeV. The parameters and are not independent; we choose V, = 1 
and then tune non-perturbatively to satisfy the lattice dispersion relation 

p^ a^p^ 
Keeping all other parameters fixed, we tune to satisfy the above relation to within ~ 1%. 

2.2 Interpolating operators for exotic and higher spin mesons 

Following Ref. [^, we employ a set of operators for exotic and higher spin mesons using 
covariant derivatives on the quark fields. In Ref. [9], such operators have forward directional 
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derivatives, V, i.e., acting only on the quark field. For states at rest, this can create required quan- 
tum numbers. However, for non-zero momentum, the forward directional derivative alone fails to 
yield operators of definite charge conjugation. One needs a backward-forward derivative V for 
this purpose. Hence we change the operator table in Ref. [Q] by replacing V — > V = V — V . 
Backward-forward derivatives were also utilized to study exotics in Ref. [S]. 



We use gauge-invariant quark smearing [g], along with stout-link smearing [IC], to obtain 
better overlap of the operators to the ground states. Smeared-source point-sink (SP), as well as 
smeared-source smeared-sink (SS), correlators are used to obtain better signal. Smearing param- 
eters are tuned separately for the various channels. We observe that displacement lengths of one 
lattice site and two lattice sites yield largely indistinguishable results in effective masses of these 
operators, and therefore for this calculation we choose displacement length one. 

3. Results 

We employ 1996 quenched configurations on a 12^ x 48 lattice (with inverse temporal lat- 
tice spacing 6.05 GeV^\ obtained from the static quark-antiquark potential. To generate meson 
correlators we used Dirichlet boundary condition with source at 5. 

In Figs. [I] - ^ we present representative effective masses for the four different combinations of 
PC, for various lattice irreducible representations Ai (0), A2(3), E{2), T\ (1) and 72(2); the numbers 
in brackets denote the lowest continuum spin lying in each irrep.. To facilitate a comparison of the 
quality of effective masses, we employ the same scale in each figure, except for the case A\ (0 ) 
where the signal is noisy and the mass larger. In each channel, including exotics but with the 
exception of the , we obtain a very good signal for several operators. Any of these operators 
would yield a good ground-state mass with reasonable statistical errors. The quality of the data is 
better than that of Ref. [^, since we adopted tuned smeared operator in part to reduce noise. A 
detailed analysis of this data using the variational method, employing multi-exponential fits, and 
using Bayesian statistics to yield ground and excited-state masses will be reported in forthcoming 
publications. 

Even our preliminary analysis enables us to make two observations. Firstly, the mass of the 
exotic 1 ^ observed in this study is lower than that found in Ref. |^, where the effective mass is 
poorer ^ We also observe a good signal for the other exotics, 0+" and 2^ . Secondly, we observe 
in the Tj^^ channel, the bottom left in Fig. 2, that one operator in this channel produces a lower 
effective mass than the others. We believe this is an overlap to the 3^ state, as T2 operators 
have an overlap to both the 2^ and 3^ states ^. In the bottom right figure in Fig. 4, we plot 
effective masses from the and A2 correlators simultaneously, and observe that they coincide 
within statistics. Finally, Fig. 5 shows the relative position of the various channels, together with 
the observed experimental ordering. 



4. Conclusions 

We present here our preliminary results on exotic and exotic state mesons in the charmonium 

'although the operators (p x B)t^ and (qq x can, at finite a, overlap with 4 ^ (a non-exotic), as well as with 
1 ^ , in the continuum no such 4 ^ can remain. As such, we are confident in assigning a 1 ^ status to the observed 
ground state. 

^An operator of the type (fcj x D)^, retains an overlap with 3^ in the continuum. 
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Figure 1: Effective masses for J^^ states with various local and displaced operators. 
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Figure 2: Effective masses for 7+ states. 



sector. We observe excellent signals for most channels, including exotics, and a comprehensive 
analysis with results for ground and excited states is in preparation. Our preliminary observation 
suggests the presence of the 1^+ exotic at around 4.2 GeV, lower than that was observed in Ref. |^]. 
The methodology developed here will be utilized in the light quark sector with dynamical clover 
fermions. Prediction of masses and particularly photo-couplings in the light quark sector will be 
quite helpful to the GlueX experiments at JLab. 
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Figure 3: Effective masses for / states. 
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Figure 4: Effective masses for J + states. Bottom right figure is a comparison of 3+ obtained from A2 and 
T2 representations. 
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Figure 5: Effective masses for various channels showing their relative positions. 
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